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Abstract 

In this paper we build models for short-term, mean-term and long-term dynamics of dune and 
megariple morphodynamics. They are models that are degenerated parabolic equations which 
are, moreover, singularly perturbed. We, then give an existence and uniqueness result for the 
short-term and mean-term models. This result is based on a time-space periodic solution exis- 
tence result for degenerated parabolic equation that we set out. Finally the short-term model 
is homogenized. 
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1 Introduction and results 

Dune and megaripple generation and dynamics, on the seabed over a continental shelf, are the results 
of interaction between the seabed and water currents. The study of the physical processes allowing 
for the generation of dunes, or governing their evolution or stability involves modeling and numerical 
simulation. Roughly speaking, the models in use essentially couple an equation for the fluid fields 
(Navier-Stokes or shallow water equations) to an equation describing sand transport on the seabed. 
Those methods were used with success in DeVriend [10J . Engelund and Hansen [11], Kennedy [19] . 
Blondeau [6 , Dawson, Johns and Soulsby [S], Johns, Soulsby and Chesher [TB], Idier Qj5] and Idier, 
Astruc and Hulsher [T7] . 

A careful watch reveals that the use of numerical simulation for the understanding of dune dynamics 
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within tide-influenced environment is essentially not efficient. The reason why is that tide oscillation 
generally prompts a coming and going of large sand volumes having a very weak resulting effect on 
dune evolution. As a consequence, questions concerning dune morphodynamics or stability have to 
be considered over large periods of time, making the computation cost expensive. 
Since many dune fields are present in strong tide region (English Channel, Celtic Sea, Irish Sea, 
North Sea, etc.) the setting out of methods to tackle dune morphodynamics in tide influenced 
environments is an important challenge. The aim of this paper is to carry out modeling methods 
and asymptotical methods for this. More precisely, we focus on linear models for seabed evolution 
and on methods which allow the removal of the explicit presence of the tide oscillations from them. 



As will be seen in section 2, for a small parameter e and constants a, b and c, equation 



~dt 



V • ((1 - &em)<? a (|u|)Vz e ) = -V • ( (1 - 6em) ffc (|u 



(1.1) 



is a relevant model for the short-term dynamics of dunes. In equation (11. ip . z e = z e (x,t) where, for 
a given constant T, t £ [0, T), stands for the dimensionless time and x — (x\, X2) € T 2 , T 2 being the 
two dimensional torus R 2 /Z 2 , is the dimensionless position variable, is the dimensionless seabed 
altitude at time t and in position x. Operators V and V- refer to gradient and divergence. Functions 
g a and g c are regular on R + and satisfy 



9a > 9c > 0, <? c (0) = g' c (Q) = 0, 
3d > 0,sup ueR+ \g a (u)\ +sup ueR+ \g' a {u)\ < d, 
sup ueR + \g c {u)\ +sup ueR+ \g' c {u)\ < d, 

3U t hr > 0, 3G t hr > 0, SUCh that U > U t hr 9a{u) > G t hr 

Fields u and m are dimensionless water velocity and height. They are given by 



where 



u(£, x) = U(t, -,x) m(t, x) = M(t, -,x), 



U =U{t, 9, x) and M = A4(t, 8, x) are regular functions on 



x T 



9 1 — > (U,Ai) is periodic of period 1, 



1^1, \M\, |^|, |^|, \VM\ are bounded by d, 

V(t,6,x) el+xRxf, \U(t,6,x)\ < U thr => 

=0, ^ = 0, VM(t,9,x) = and VU(t,d,x) = 



dt ' dt 
39 a <9 u] e [0, 1] such that V 9 e [0 a ,9 u ] 



\U(t,9,x)\ > Uthr- 



(1.2) 



(1.3) 



(1.4) 



Remark 1 . 1 The last two assumptions in (11.41) are necessary when g a may vanish. In the case where 
9a{u) > G t hr for any u > 0, then U t hr = and the last two assumptions of (|1.4I) are automatically 
satisfied by any U. ■ 



The following equation, for constants a, b and c 

r)z^ CI CI 11 \ 

_ _ _ V • ((1 - 6^m) 5Q (|u|)Vz e ) = -V • f(l - 6^m) 5c (|u|) — J , (1.5) 
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with condition (|1.2[) on g a and g c and with u and m given by 

u(t,x) = U(t ) ^=,-,x), m(t,x) =M(t,-=,-,x), 
ye e e 

is a relevant model for mean-term dune dynamics. 
For mathematical reasons, we assume 



(1.6) 



U(t,T,9,x) = U(t,9,x) + ^Ui{t,T,9,x), 



(1.7) 



where U = U(t,9,x) and U\ — Ui(t,T,0,x) are regular. We also assume that M — M(t,r,9,x) is 
regular and 



i — > (U,Ui,M) is periodic of period 1, 
r i — > (Ui,M) is periodic of period 1, 

PI l-^l. 1^1, |VW|, |Wi|, |-^|, |-^|, |VWi|, 

S.M 9.M dM 

\M\, | — |, | — |, | — |, | VX | are bounded by d, 

V(t, t, 9, x) G M + xRxIx T 2 , 0, i)| < f7 t/lt . =^> 

— (t,e,x) = o, —(t,9,x) = o, vu(t,T,e,x) = o, 

dt dM dT dM 

—-(t,9,x) = 0, -z—(t,9,x) = and VM{t,6,x) = 0, 
at ot 

36> Q < 9 U e [0, 1] such that V 6 e [6 a> 6 u ] => \U{t, r, 0, x)\ > U thr . 



(1.8) 



A relevant model for long-term dune dynamics is the following equation 



^ _ £ V ((1 - bem)g a (\u\)Vz*) = • f (1 - bem)g c (\n\)^ 



(1.9) 



where a, & and c are constants, where g a and g c satisfy assumption (|1.2[) , and where z is defined on 
the same space as before. It is also relevant to assume 



u(x,t) = U(-,x) + t 2 U 2 {t, -,x), m(t,x) = M(-,x) + e 2 M 2 (t, -,x) 



(1.10) 



where U = U{0, x), U 2 (t,9,x), M = M(9, x) and M 2 = M 2 (t, 6, x) are regular and 



9 i — > (U,U 2 ,M,M 2 ) is periodic of period 1, 
dU dU dU 2 ,dU 2 dM. dM. .dM. 



<9t " 1 96> 



dr 



dt "'m " -» '' TJ at 

l^al, l^r 1 !, |^|, 1^1, |VX 2 | arc bounded by d, 



dt " 1 90 " 1 dr 
V(t,9,x) el+xMxT 2 , \U(t,9,x) + e 2 U 2 {t,8,x)\ < U thr =>■ 

— -(d,x) = 0, VW(0,s)=O, VW 2 (0,:r) = 0, 
dt dM, 

= 0, VM(0,x) = 0, VM 2 (i,S,i) = 0, 
30 a < W e [0, 1] such that VfleM, e [0 a , W ] =*■ \U(9, x) + e 2 W 2 (i, 9,x)\> U thr 



(1.11) 
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Equations (jl.lj) . (| 1 . 5|) or (|1.9I) need to be provided with an initial condition 

zf t =o = *b. (1-12) 
giving the shape of the seabed at the initial time. 

We will now state the main results of the paper. The first concerns existence and uniqueness for the 
short and mean-term models. 

Theorem 1.1 For any T > 0, any a > 0, any real constants b and c and any e > 0, under 
assumptions U.ty) and \1.S\) , or \l.T\j, \1.8\) , if 

z Q e L 2 (T 2 ), (1.13) 

there exists a unique function z e G L°°([0, T), L 2 (T 2 )), solution to equation M.l)) or hl.5\) provided 
with initial condition U.12\ ). 
Moreover, for any t 6 [0,T], z e satisfy 

lk C ||L = ([0,T),L 2 (T 2 )) < 7) (l-l 4 ) 

/or a constant 7 wot depending on e and 



d ( / z e (i, x) dx 

/T 2 

dt 



= 0. (1.15) 



The proof of this theorem is done in section [3] except equality (j!.15[) which is directly gotten by 
integrating ([l.ip or (|1.5p with respect to x over T 2 . 

In the previous theorem, L 2 (T 2 ) stands for the usual space of square integrable functions defined on 
the torus T 2 and L°°([0, T), L 2 (T 2 )) stands for the space of functions mapping [0,T] to L 2 (T 2 ) and 
which are bounded. ||-||l°°([o.t),l 2 (t 2 )) stands for the usual norm on this space. 

Remark 1.2 As equations (ll.lj) and (jl.5[) are linear, almost parabolic equations, the proof of the 
existence of z e over a time interval depending on e is a straight forward consequence of adaptations 
of results from Ladyzenskaja, Solonnikov and Ural'Ceva [2T] or Lions [22]. But here, since we want 
to follow an asymptotic process consisting in making e — > 0, we need a time interval which does not 
depend on e. Because of the presence of - factor and the fact that the diffusion term may cancel, the 
proof of theorem 11.11 needs several steps. In a first step, we prove the existence of a solution, periodic 
in time and space of a parabolic equation. From this first existence result, we deduce existence of a 
solution, periodic in time and space of an ad-doc degenerate parabolic equation. 
Those two results are interesting by themselves and complete the theorem collection in the topic 
of time and space time-periodic solution to parabolic equation. Concerning this topic, we refer for 
instance to Barles and Souganidis [3J, Berestycki, Hamel and Roques [4j [5] , Bostan [7], Hansbo |15) . 
Kono [20], Nadin [25] [24], Namah and Roquejoffre [26] and Pardoux [29] . 

Then, having on hand the existence of the space-time periodic solution to the ad-doc degenerate 
parabolic equation, we can deduce that the solution z e which exists on e-dependant time interval, 
remains close to it. This allows us to deduce a large time existence. ■ 

Remark 1.3 Moreover, notice that theorem theorems 13.161 and 13.171 also complete the theo- 
rem collection concerning the topic of large time behavior of parabolic equation (see Barles and 
Souganidis [3 , Da Lio [8], Norris [28], Park and Tanabe [30], Pardoux [29], Petita [31] and Tanabe 
[32]. ■ 
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We now give a result concerning the asymptotic behavior as e — > of the short-term model. 

Theorem 1.2 For any T > 0, under the same assumptions as in theorem \l.ll the solution z e to 
equation i f 1.1]) given by theorem ] 1.1\ two-scale converges to a profile U G L°°([0, T], X^(R, i 2 (T 2 ))) 
which is the unique solution to 

^L-V.(AVU)=V-C, (1.16) 

where A and C are given by 

A = ag a {\U{t,6,x)\) andC = cg c (\U(t,6,x)\) ,^'1r ( L17 ) 

\U(t,8,x)\ 

In this theorem, L2'(M, L 2 (T 2 )) stands for the space of functions depending on # and x mapping 
R to L 2 (T 2 ) and which are periodic of period 1 with respect to 6 and £°°([0, T],L^(R, i 2 (T 2 ))) 
stands for the space of functions mapping [0, T] to L^(R, L 2 (T 2 )) and which are bounded. For the 
definition and results about two-scale convergence we refer to Nguetseng [37] , Allaire [T] and Frenod 
Raviart and Sonnendriicker |13j . ■ 

Finally, we give a corrector result for the short-term model under restrictive assumptions. 

Theorem 1.3 Under the same assumptions as in theorem \l.l\ and if moreover Uthr = 0, considering 
function z € € L°°([0, T), L 2 (T 2 )), solution to U.l\) with initial condition \1.12\) and function C/ £ g 
X°°([0,T],L^(]R,L 2 (T 2 ))) defined by 



U e (t,x) = U(t, (1.18) 



where U is the solution to (1,16\) . the following estimate is satisfied: 

z e - V 



<a, (1.19) 

L~([0,T),L 2 (T 2 )) 



where a is a constant not depending on e 
Furthermore 

z e - V 



two-scale converges to a profile U 1 G L°°([0, T],L°»(R, L 2 (T 2 ))), (1.20) 



r 

which is the unique solution to 

dU 1 „ /t_„i\ _ -x dU 



V • [AVU 1 J = V • Ci + — + V • (A\VU), (1.21) 



de v j dt 

where A and C are given by j!.17\) and where A\ and C\ are given by 
Ai(t,9,x) = -abM(t,6,x)g a (\U(t,6,x)\), C x {t,0,x) = -cbM(t,6,x) g c (\U(t,6,x)\) U ^ 6 > x ) 



\U{t,9,x)\ 

(1.22) 



Remark 1.4 Theorems 11.21 and 11.31 state a rigorous version of asymptotic expansion of z e : 

z € (t,x) = U{t, t -,x) + eU 1 (t, t -,x) + ... . (1.23) 

■ 
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2 Modeling 

2.1 Sand transport equation 

The equation modeling sand transport is the following (see Van Rijn [33], Idier [16]): 

dz 1 

g+_^V-«z = 0. (2.1) 
at 1 — p 

In this equation the fields depends on time t € [0, T), for T > 0, on the horizontal position a; = 
(xi,X2) G f2, where is a regular open set of K 2 . The field z = z(t,x) is the height of the seabed 
in position x and at time t and g = g(a;, t) is the sand volume flow in x and at t. The parameter 
p G [0, 1) is called sand porosity. Equation (|2.ip has to be coupled with a low linking the sand flow 
q with the seabed height variation and the velocity of the water near the seabed. Usually, such a 
law is written 

q = q f -\q f \XVz, (2.2) 

where g/ stands for the water velocity induced sand flow on a flat seabed and where |g/| stands for 
its norm. The constant A is the inverse value of the maximum slope of the sediment surface when 
the water velocity is 0. A generic way to write qj is 

If = ax(s(H) - ff(0)A , (2.3) 

where g is a non-negative regular function defined on K + and where x is a regular function from R 
to R, being on R~ and increasing on R + . u is the water velocity near the seabed, g(u) is regular 
function of u £ ]R + and u c is the threshold under which the water velocity does not make the sand 
move. 

Every law encountered in the literature, for instance Meyer-Peter and Miiller |23] formula, Bagnold 
and Gadd formula (see [2] and [14]) and Van Rijn [33] formula, is recovered by setting functions \ 
and g. 

In the sequel of the present paper we shall restrict ourselves to laws of the Van Rijn type [33 which 
consists in writting 

q f =a X (D G (\T b \-T c ))^-, (2.4) 
where r b is the shear stress density imposed by the water on the seabed. It is linked with u by 

r b = P^ R , (2.5) 

where p is the water density, C is a constant defined by C = l n (^5^-), d being the water height 
above the seabed and Dq being the sand speck diameter. The threshold r c expresses as 

T C = P§ I , (2.6) 

and x is given by 

y(<j) =0 if a < 0, 

V ' (2.7) 

= |(7 3 / 2 | if(T>0. 

The order of magnitude of constant a is 100. 
Injecting equation (|2.5p into (|2.4j) and (|2.2j) we get 



u c 2 \ I u 



q = a X [D GP ' 1 c2 ) ( — -AVz) , (2.8) 



u 
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and equation (12.11) reads 

dz a 



dt 1-p 



V • 



X\ D GP- 



c 2 



= o. 



(2.9) 



2.2 Scaling 



Now, we will scale (|2.9[) to write a dimensionless version of it. We introduce a characteristic time t 
and a characteristic length L and we define the dimensionless variables t' and x', making t and L 
the units by 

t = it', x = Lx'. (2.10) 



We also define z the characteristic height of the dunes and the dimensionless seabed height 

z'(t' ,x') — —z(tt',Lx'). 



(2.H) 



Concerning coefficients of equation (|2.9p . we introduce u the characteristic velocity of the water, we 
consider the mean water height H and M the characteristic height variation due to the tide. Then 
we define u' being the dimensionless water height variation by 



u'(i', x 1 ) = -u(tt', Lx'), m'(t', x') = -=-{d(tt'Lx') - H). 
u M 



(2.12) 



Once those variables and fields are introduced, we first approximate C, taking into account that ^ 
is small. 

c = "<^) +1 °( 1 + f m ' |; 



From (|2~T3l) we get 



1 

C3 



hi 



AH 



1-3- 



f4H\ M , 



(2.13) 



M 



AH 



Since for instance 



Vz{tt',Lx') = —Vz'(t',x'), 
zL 



(2.14) 



(2.15) 



we get from equation (|2.9p the following equation for z' 



dz!_ 

dt r 1-p 



M 



1 iu 3 ( P D G )V 2 



.'|2 



V'z' 



1-p 



1 - 3- 



M 



u 



(H^)) c Lz 

. (2 - 16) 

Having this dimensionless model on hand, we will now consider several situations in setting the 
characteristic values for short, mean and long-term dune evolution and for small and big sand 
specks. 



First, we fix the characteristic sizes which are common for every situation. Dunes exist within coastal 
ocean waters over a relatively flat continental shelf, with a water height of about 30 to 50 meters, 
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with tide induced height variations which are not too strong and with relatively strong tide currents. 
Then we set 

u = lm/s, H = 50m, M = 5m. (2.17) 
Moreover, the order of magnitude of coefficient is 1, then we get 

= 1 and — !— = 2. (2.18) 
1-p 1-p v ' 

Now we detail the sizes of every characteristic value and of their concerned ratios in equation (|2.14l) , 
(|2.16|) for every situation. 

Short-term dynamics of dunes made of a small sand specks 

Here, we shall consider that i is an observation period of time. We take as t the order of magnitude 
of the smallest period of time during which dunes undergo significant evolution in a tide-submitted 
environment, i.e. i = 100 days ~ 2400 hours ~ 8.6 10 6 s. Introducing ui the main tide frequency, t 
has to be compared with the main tide period A ~ 13 hours ~ 4.7 10 4 s. This leads to the definition 
of a small parameter e: 

— = e. (2.19) 
tu 200 v ; 

We consider that the sand speck diameter Dq is 0.1mm = 10~ 4 m. According to Flemming [T2] and 
Idier [16], this gives rise to dunes being about 1 meter high, the wave length of which is about 10 
meters. Then we set 

z = lm and L = 10m. (2.20) 
We also consider that the critical velocity u c is small compared with u. In other words we set 

2 =0. (2.21) 

u A 

As the computations of the factors in (I2.20[) yields 



1-P f^(4H^ 3 l2 2e' 



A „ iu 3 (pD G ) 3/2 . 1onn . 10 



-a " r 1800 , (2.22) 

1-P (^f^lz 



3M 



2.10" 2 - 4e, 



equation (|2.16l) reads 

|£ _ I V • ((1 - 4em)|u| 3 Vz) - • ((1 - 4em)|u| 2 u), (2.23) 
where we removed the '. 

Concerning fluid fields u and m, we assume that they are periodic functions, with modulated am- 
plitude, and of period the tide period. In other words 

n{x,t)=U{t,-,x), m(x,t) = M{t,-,x), (2.24) 

for functions U and A4 being regular, and such that 8 i — > (W(£, 9, x),A4(t, 9, x)) is periodic of period 
1, with a null mean value. 
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Finally, as dunes of the considered kind are, in nature, gathered into dunes fields it is not completely 
unrealistic to set equation (j2 .23[) in a periodic position space. 

As matter of the fact, considering equation (|2.23p is appropriate for the study of short-term dynamics 
of dunes made of small sand specks with a mathematical point of view. 



Short-term dynamics of dunes made of a big sand specks 

For this regime, we consider: 

t ~ 100 days ~ 2400 hours ~ 8.6 10 6 s, 



- ~ 13 hours ~ 4.7 10 4 s Do = 5 10~ 3 m (2.25) 

W - 1 

z — 50 77i, L = 300 m, u c — — m/s. 



Then 

A ftW/* 1 



A iu 3 ( P D G ) 3 / 2 5 , x 
" 3 1000 --, (2.26) 



1-P (^^L<\- L -z 



3M 



1.3 10" 2 - 3e, 



•Do 

equation (|2 . 1 6[) . with the ' removed, gives 



g - lv • ((1 - 3em) X (|u| 2 - i)V*) - Jv ■ [(1 - 3 £ m) X (|u| 2 - 5) AY (2-27) 



Mean-term dynamics of dunes made of a small sand specks 

By mean-term we mean a period of time of 4.5 years ~ 54months — 1.4 10 8 s. Then, we take 
t = 1.4 10 8 s, which is compared with -1 — 13 hours — 4.7 10 4 s giving 

k ~ 30CX) - ' (2 - 28) 

We also consider a second tide period which is the time for the earth, the moon and the sun to 
recover approximately the same relative positions. This period of time J- is about one month. So 
we have 

7^ ~ ~ VI (2-29) 
tuj c 54 

We also take Dq — 5 10~ 5 m and 

z = l 777, L = 10 777, u c = 0. (2.30) 
Computing the coefficients in equation (|2.16[) gives 

ft? 1 20 

°JL _ i V • ((1 - V^m)|u| 3 Vz) = -V • ((1 - ^m)|u| 2 u). (2.31) 

at e e 

As was previously seen, it is reasonable to set this equation in a periodic domain and concerning 
the fluid fields we consider 

u(t, x) = U(t, -p, -), m(t, x) = M(t,-^=,-). (2.32) 
V e e v e e 
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to take into account the two tide periods under consideration. In (|2.32p we take U and M. as regular 
functions such that 

r i — > (U(t, t, 0, x),M(t, t, 9, x)) 



(2.33) 

(U{t,r,e,x),M(t,r,e,x)) 



are periodic of period 1. 



Long-term dynamics of dunes made of small sand specks 

We take here f ~ 16 years ~ 1.4 10 5 hours ~ 5 10 9 s. We compare this period of time with the 
second tide period ~ 1 month ~ 2.6 10 6 s. Then, we define e by 

1 1 

== = e. (2.34) 

tu c 192 v ' 

We set 

D G = 7.10~ 5 m, z = lm, L = 10m, u c = m/s. (2.35) 
with those values equation (|2 . 16[) yields 

dz 1 20 

— + T V • ((1 - 4em)|u| 3 Vz) = -5-V • ((1 - 4em)|u| 2 u). (2.36) 
ot e z e z 

As, at the second tide period scale the tide phenomena may almost be considered as really periodic 
we set 

u(x,t) = U{-,x) + e 2 U 2 (t, -,x), 

t % ( 2 ' 37 ) 
m(a;,t) = M(-,x) + e 2 M 2 (t, -,x), 
e e 

where U, U2, M M2 are regular functions such that 8 1 — > (U(0,x),U2(t,0,x),M.(9,x),M.2(t,9,x)) 
is periodic of period 1 and such that 

1 

U(9,x)d9 = 0, (2.38) 


1 

M(6,x)d6 = 0. (2.39) 

3 Existence and estimates, proof of theorem 11.11 

Defining 

A £ {t,x) = Ae(t,-^=,-,x), (3.1) 

where 

A e (t,T,6,x) = a(l - bs/eM(t,T,0,x)) g a (\U(t,6,x) + ^/eUi(t,T,9,x)\j, (3.2) 

C e (t : x) =C e {t,-^=,-,x), (3.3) 

and where 

C e (t,r,9,x) = c(l - b^M(t 7 T 7 9,x)^ g c (\U(t,e,x)+y/eUi(t,T,6,x)\' 

U(t,8,x) + y/eU(t, t,6,x) 



\U(t,6,x)+«ftU(t, r,0,z)r 



(3.4) 
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equation (|1.5p . (|1.12[) with assumptions (|1.6[) and (|1.7p reads 



^---V-(^Vz e ) = -V-C e , 
ot e e 

Z\ t=0 = z Q . 
In the same way, setting 

A e (t,T,9,x) = A e {t,9,x) = o(l - &e.M (t , 0, a;) ) g a ( \U (t , 9, x) | ) , 

and 



C e (t, t, 9, x) = C e (t, 9, x) = c(l - beM(t, 9, x)) g c (\U(t, 9,x)\) 



U{t,9,x) 

\u(t,e,x)\' 



(3.5) 



(3.6) 



(3.7) 



and defining A e and C e from A c and C e by (I3.1[) and (I3.3p . we may deduce that equation (jl.ll) . with 
assumption (jl.3l) . can be set in the form (|3.5p . 



From assumptions (fTT2|l and pT4|) or (|L2j) and (|L8j) . „4 e defined by (pT2|t or (|3J|) and C e defined by 
p.4p or (|3.7p satisfy the following properties 



|A| <7, |C e |<7, 



dA e 



<7, 



90 



dr 



at 



dA t 



dt 



<7> 



dC e 



Of 



<7, 



<7, |VA|<7, |V-C £ |<7, 
dV-C f 



< 



ac. 



dt 



< 



7: 



dr 



< Ve7) 



3r 



< ^7, 



(3.8) 



(3.9) 



on R + x t x 1 x T 2 , for a constant 7 depending only on a, 6, c and d and not on e. 
Concerning (13.91) in the case where A e and C e are defined by (|3.6I) and (13.71) , it reduces to 



8A e dA e 



dr di 



= 0. 



(3.10) 



Moreover, for every e, < e < 1, A\ > 0, 



r 1 — ► (_4 e ,C e ) is periodic of period 1, 
1 — ► (,A e ,C e ) is periodic of period 1, 



(3.11) 



and there exists a constant Gthr depending only on a, b, d and Gthr and two numbers 9 a and 0^ in 
[0, 1], 9 a < 9 U , such that 

A t (t,T,9,x)>G thr , (3.12) 
for every t G R, r G R, x G T 2 and 6> G [0 a , W ] and such that V(i, T,»,i)el+xlxlxT 2 



A e (t,r,9,x) <G 



Ihr 



dA 8 A ~ 

-£-(t,T,e,x) = Q, -j^(t,r,6,x) = 0, VA(i,r,M) = 0, 
at dr 

dC dC ~ 

-^(t, T , 9, x) = 0, -^{t, T , e,x) = o,v- e £ (t, t, 9, x) = 0. 



(3.13) 
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We also have the following inequalities 



\Ce\ <l\A e \, 
|Ce| 2 < 7 |A| 

|VA| < 7 |A 
dA 



dt 

S(VA) 2 



dt 
8A ( 2 



dr 

OVAr 2 



dr 
0C f 



<7|A|, 

<7|A|, 
<£7|A|, 
<£7|X|, 



9t 

2 



<7|A|, 
<7|A|, 

<7 2 |A|, 



3.14) 
3.15) 
3.16) 

3.17) 

3.18) 

3.19) 

3.20) 
3.21) 

3.22) 

3.23) 



possibly changing the value of 7 and making it also depend on Gthr ■ 

Inequality (I3.14|) is a direct consequence of (|1.2[) . Since for small values of |C £ |, |C e | 2 < |C £ | and since 
C e and A e are bounded, inequality (|3.15[) follows from p,14p . When |y^ e | < G t hr, then \7A e = 0. 
Hence (|3.16[) is realized. When |_A e | > G t hr, since ^4 £ and \7A e are bounded, p,16p is obviously 
realized. Hence f|3 . is true. With a similar argument (|3.17l) - (|3.20l) may be obtained. In order to 
obtain (|3.21[) . we just have to notice that, when \A e \ < Gthr, V • C e = 0. Hence we can give the same 
argument as above. In the same manner, the last two inequalities may be obtained. 



We now consider for a positive small parameter v the following regularization of p. 5 



Q z e,u x 



dt 



V • {{A' + v)Vz e > 



1. 



-v-c e 



(3.24) 



-|t=o 



z . 



Denoting by || • H2 and || • Hoc the usual norms of spaces L 2 (T 2 ) and L°°(T 2 ), applying the energy 
estimate and the maximum principle, (see for instance Lazyzenskaja, Solonnikov and Ural'Ceva|21] 
or Lions [32]) we can obtain the following lemma. 

Lemma 3.1 For any T > 0, if z e L 2 n L°°(T 2 ) under assumptions H3.ll}) , then for any e > and 
v > 0, there exists a unique solution z e ' v £ L°°([0, T); L 2 n i°°(T 2 )) to [3.2$ . Moreover it satisfies 

■ € '"\\oo < -, (3.25) 



l^lla 



for a constant 71 depending only on 7 and \\z0W2 + H^olloo- 

As estimate (|3.25p depends on v, letting v go towards we obtain the following corollary. 
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COROLLARY 3.2 For any T > 0, if z e L 2 n L°°(T 2 ), and under assumptions h3.ll}) , then for any 
e > 0, there exists a unique solution z e £ L°°([0, T); L 2 n L°°(T 2 )) to (3.5\) . Moreover it satisfies 

lk e l|2 + lkloo< ( 3 - 26 ) 

The uniqueness of the solution to (13.5[) is a direct consequence of the linearity of this equation. 



As we want to study the asymptotic behavior of z e as e goes to 0, estimates (|3.26p and (|3.25[> arc 
not enough. We need estimates which do not depend on e. For this, we first consider the following 



problems, which consists in finding S v 
1 with respect to 9, solutions to 



S"(t, r, 6, x) and = S"(t, r, 6, x) being periodic of period 



— - V • ((A e (t, r, •, •) + i/)V5") - V ■ C e (t, r, •, •), 



and 



V • r, •, •) + ^)V<0 = V • C e (t, r, •, •)• 



(3.27) 



(3.28) 



In equations (13.271) and (|3 . 28[) t and r are only parameters. 

The method to get the desired estimates which do not depend on e is shared in several steps. In the 
first, we set out the existence of periodic solution to equation (I3.28|) . We, moreover, set out that 
sequence is bounded independently of fi and e. We also show that 5^ is differentiable with respect 
to t and t. In a second step, letting u, go to 0, we get existence of S u , with the same properties as 



S". The third step consists in finding estimates on S v 



US' 
dt 



and $f- 

OT 



which are independent of v to 



be able to make the process v — > and to obtain the existence of a solution S to (|3.27|) with v = 
and consequently a periodic solution Z e to an equation close to (|3.5[) in a fourth step. The fifth step 
consists in noticing that the solution z e of (|3.5|) is not far from Z e . 

The framework of periodic solutions to parabolic equations is widely studied in both linear and 
nonlinear cases. We refer for instance to Barles[3], Berestycki-Hamel and Roques [5j|4], Bostan [7], 
Hansbo [15], Kono [20], Nadin [25j[24], Namah and Roquejoffre 26, and Pardoux 29! for a revue of 
the results on this topic that our result (see theorem 13 . 3 1 and 1 3 . 1 j) completes. Inspired by ideas that 
may be found in those references, concerning equation (|3.28|) we can state the following theorem. 



Theorem 3.3 Under assumptions 113. 8\) . (3. 9\) and 113.11)) . for any /i > and any v > 0, there 
exists a unique = S"(t,T,0,x) G C° l~l L 2 (M. x T 2 ), periodic of period 1 with respect to 8, solution 



to S3.28\) and regular with respect to the parameters t and r. Moreover, there exists a constant 73, 
which depends only on 7 and v such that 



sup 

8ER 



It 2 



S u J9,x)dx 



0, 



de 



< 73, 



liv^iu- 



B,L 2 (T 2 ) < 73, 
!,L 2 (T 2 ) < 73- 

The following estimates with respect to the parameters t and r are also true 



dt 



< 73, 



1,L 2 (T 2 )) 



dr 



L~(R,L 2 (T 2 )) 



< Ve73 



(3.29) 

(3.30) 

(3.31) 
(3.32) 

(3.33) 
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In the above theorem, the norms are defined by 



ll/llk (RXT3) = / / fdxdB, (3.34) 

* Jo JT 2 

ll/lli~(R,L2( T 2)) = sup / fdx. (3.35) 
* ee[o,i]JT 2 

Proof of theorem 13.31 The point of departure to prove theorem 13.31 consists in considering, for 
£ G L 2 (T 2 ), the solution to 

< + §-V-((A + ,)VQ = V. C ~, (3 36) 

whose existence and uniqueness on any finite interval is a direct consequence of Ladyzenskaja, 
Solonnikov and Ural'Ceva [2T] or Lions [22]. We also consider the application 

□ : L 2 (T 2 ) —> L 2 (T 2 ), £^^(1,.), ( 3 - 37 ) 

and for it we can prove the following lemma. 

Lemma 3.4 For any \i > 0, v > and e > 0, □ is a strict contraction 

Proof of lemma [3.41 For any £ 6 L 2 (T 2 ) and any £ S i 2 (T 2 ), we consider £^ and £^ the solutions 
of (|3.36[) associated with initial conditions £ and £. It is obvious to obtain that £^ — £^ is solution to 

M(C - O + g 9 - V ' ((A + "M^ ~ O) = °> ( 3 - 38 ) 



which multiplied by £^ — £^ and integrated on T 2 yields 



/'llC-Clll+a-flfl + JJ^ + ^mZ-O \ 2 dx = 0. (3.39) 

Since A t + v > 0, we can deduce that 

11^(1, ") - ^(1, < e- 2 ^||C - elll, (3-40) 
giving the lemma. ■ 

From lemma l3~4l we deduce that there exists a unique function £ 6 L 2 (T 2 ) such that □(£) = C- 
Since A t and C e are periodic with respect to 9 we deduce that the sought periodic solution of 
(|3.28[) is nothing but the solution of (|3.36[) associated with initial condition £ such that □(£) = C- 
Hence we proved the existence of claimed in theorem 13. 31 The fact is that is continuous comes 
from the fact that it is a solution of a parabolic equation with regular coefficients. We now turn to 
the properties of the function S". 

Lemma 3.5 For any u, > 0, v > and e > the solution to S3. 28\) satisfies property L3.29\) . 
Proof of lemma 13.51 Integrating (|3.28|) over T 2 gives 

dff s;dx 
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which leads, 



T 2 



S"(6,x)dx = / S v J6,x)dx e-^ 



J 2 



(3.42) 



Since is periodic of period 1 with respect to 8, J T2 S^dx is also periodic of period 1. Hence the 
only possibility to satisfy (13.42[) for J T2 S^dx is to be 0. This ends the proof of lemma 13751 ■ 



Lemma 3.6 For any /i > 0, v > and e > estimate 

l|V<S^|| L |( KxT 2) < -, 



and estimate II 3. 3 Op are valid. 

Proof of lemma 13.61 Multiplying equation (|3.28[) by S" and integrating on T 2 give 



tA\s v M-)\\i 



2 , ljOTOOjlj) 



T- 



d0 jf2 

v -c e s»(e,-)dx 



(A e + v)\VS»(6,-)\ 2 dx 



e 6 vs;(9,-)dx. 



Integrating (|3.44p . from to 1, with respect to 9 gives 

-l 



v\\2 

AillL^(RxT 2 ) 



JT 



{Ae + v)\VS^\ 2 dxdO 



10 JT 2 

Since A e + v > v, we deduce from (|3.45|) 

-.i/i|2 



C e V^(^-)^<7||V^|U |(RxT2) . 



HI V ^|| L 2 (RxT2) ^ 7||V5^|| L 2 (RxT2) , 



giving f3~43l) . 

Multiplying (13.28)) by and integrating on T 2 give 



d\\S v M 
1 dd 



Oft 



Since 



d\ J {A c + v)\VS^\ 2 dx 
d8 



~ dS v f ~ dS v 

jX + ^s;.v^dx = lv.c^dx 



^\vs;\ 2 dx + 2jjA t + ^s;.^dx, 



integrating (I3.47[) with respect to 8 from to 1 gives 
2 -i 



Oft 



^s;\ 2 dxdo 

L|(RxT2) JO JT 2 00 



JT 



BC 

2 -^VS'dxdO 



< 7(I|V5^|| L 2 (RxT 2) + \\^S^\\ L 2 #{RxT 2)) 



Because of (|3.43[) . we can deduce that (13.30[) is also true, ending the proof of lemma 



(3.43) 



(3.44) 



(3.45) 



(3.46) 



(3.47) 



(3.48) 



(3.49) 
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Lemma 3.7 For any n. > 0, v > and e > 0, the solution to 13.28\) satisfies 



(3.50) 



and estimate (3.31]) . 



Proof of lemma 13.71 Multiplying (13.281) by — A6>^ and integrating with respect to x G T 2 yields, 



Mlivs 



„.. a , ld(HV^Hl) 



r/(9 



V • ((A + i/)VS!;) ASfJda: = - / V • C 6 A5"da;, (3.51) 



or 



A*I|V5: 



/ M I 2 



ld(l|V5-|||) 



T 2 



Since for any real number U and V, 



\ UV \ < ^ll^i/i 



T 2 



1 



T 2 



4 A + f 

using this formula with U — AS^ and V = VA • V<S^ we get 



(A + u)\AS%\*dx = - I VA • VS^AS^da; - / V • C e AS^dx. 

(3.52) 

(3.53) 



J 



VA ■ WS U U AS u u dx 



< 



T 2 



A ' + ~\^dx 



|VA| 2 

T 2 A + v 



In the same way, applying Q3.53[) with {7 = A6>^ and U = V • C e we get 



(3.54) 



T 2 



V • C e AS v a dx 



< 



^^\AS^dx- 



T 2 A + ^ 



-dx. 



Using (|3T54l and (j3~55j) in (|3~52"j) gives 

M||V^|| 2 + - 

which, integrating in from to 1, yields 



A_+_^ 

T 2 2 



\As;\ 2 dx<^(( \ws;\ 2 dx + i), 

V Jj2 



^W^S'iIWlI (kxt 2 ) 



A + f 



o Jt 2 



and in particular 



ll AlS Mlli^(RxT 2 ) 



2 

T 2)< r (- + 1), 



leading to (|3~5TJ1) . 

From lemma [3~7J1 we may deduce that there exists a Qq G [0, 1] such that 

iiv<m,-)ii 2 <?. 



(3.55) 



(3.56) 



\AS»\ 2 dxdt < i-(||V^||| 2#(RxT2) + 1), . (3.57) 



(3.58) 



(3.59) 



Beside this, from f|3 . 56[) we can deduce 



d(\\VS, 



a* lla) 



27 2 



d0 



<— / ivs^ + i 



(3.60) 
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which gives, integrating from 9o to any other Q\ £ [0, 1], 

HV^((?i,.)||2-||V5^o,-)||I<^7 f 1 (I \VS^dx + l)d6 

v Je \Jt 2 ' / 



<^(l|V^||i |(RxT2) + l), (3.61) 



giving the sought bound on ||V«S„(0i, Olli f° r an Y ®\ or m other words (|3.31[) (possibility changing 
the constant 73.) ■ 



Lemma 3.8 For any u. > 0, v > arid e > £/ie solution to \3. 28\) satisfies estimate (3.32)) . 
Proof of lemma 13.81 At any 8 G R, we consider the Fourier expansion of <S^(0, -) : 

S*(0,-) = J2$ne mx . (3.62) 



ner 

We have 



\\s;(e,-)\\l= E i^i 2 - ( 3 - 63 ) 

nEN 2 

HV5^,-)lla= E N 2 I^«| 2 ' ( 3 - 64 ) 

nSN 2 

and, as a consequence of (|3.29[) 

So = 0. (3.65) 
Then, we obviously deduce from (I3.63|) - (l3.65p that 

K(^-)||I<||V^(^-)||I, (3.66) 
for any 9, giving (|3.29p as a consequence. This ends the proof of lemma [ 



Finally, we can prove the following lemma 

Lemma 3.9 For any [i > 0, v > and e > the solution to VS. 28\) satisfies estimate 



Proof of lemma I3.9L Obviously, we get that is a solution to 



d 



""af + do v ■ ( ( ^' T ' ■' + v) v ^f ) = v " 1 ■ ■ i:U)7) 



with 



(3.68) 



From the previous estimates, we can deduce that there exists a constant 75, which only depends on 
7 and v, such that 

||Ce|U|( Rx T=) <75, (3 _ 6g) 
V ' C e ||i2 (RxT 2 ) < 75- 
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Proceeding in the same way as when proving lemma [3~51 we get 



o Jt 2 



(A, + v) 



dsr, 



at 



dxdO < 75 



V — - 

dt 



L|(RxT 2 ) 



giving 



Of 



< 



75 



L|(RxT 2 ) V 



9S V U 
V — - 

dt 



V- 



dt 



7l 



d0 



< — ( 



;et 








/ 


V — - 


It 2 


dt 



Ixf 2 ) 



1), 



and using (|3.70p , from which we can deduce 

dsr, 



<2!L. 

2 1/ 



for a given 9q, we finally obtain 



V M 



0t 



.,L 2 (T 2 )) 



is bounded. 



Since the mean value of -§£-(6, •) is for any 9 6 R, using the same argument as in lemma 15781 we 
conclude that the first estimate of Q3.33P is true. We can do the same for the second estimate, to 
end the proof of the lemma and also of theorem [ 



(3.70) 



(3.71) 



(3.72) 



(3.73) 



(3.74) 



Since the estimates in theorem 13.31 do not depend on fi, making the process fi — > allows us to 
deduce the following theorem. 



Theorem 3.10 Under the assumptions Iff.fflj .i fXP)) and hS.ll}) , for any v > 0, there exists a unique 
S u = S v (t, r, 9, x) £ L 2 (lxT 2 ), periodic of period 1 with respect to 9 solution to \3.27\) and submitted 
to the constraint 



sup 



S"(9,x)dx 



= 0. 



Moreover, there exists a constant 73 which depends only on 7 and v such that 

dS v 



09 



<T 2 ) 



< 73, 



\\VS U \\l^(M,L 2 (T 2 )) < 73, ||<S1l^(r,l 2 (t 2 )) < 73, 



dS u 



dt 



< 73, 



dS" 



di 



< \/e73- 



(3.75) 

(3.76) 
(3.77) 



Proof of theorem 13.101 As previously said, existence of S v follows from making /z tend to in 
(J3HH1)- Formulas f3~75|) - (j3~77|) directly come from estimates 

Uniqueness is insured by p. 751) . once noticed that, if S v and S v are two solutions of p.27|) . with 
constraint (I3.75|) . S v — S v is solution to 



d{S u - S v ) 
d9 



V • ((A + i/)V(S" - S v )) = 0, 



(3.78) 
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from which we can deduce that 

H|v(^-^)|||^ (RxT2) = o, 

and because of (|3.75p . and its consequence: 

IIS" - <S' y ||L^(R,T 2 ) < 11^(5* ~ S Iy )|| i 2 t(RxT 2 ) , 



that 



(3.79) 

(3.80) 
(3.81) 



This ends the proof of theorem 13.101 



Having on hand theorem 13.101 we will set out the properties of 5" which will allow us to make the 
process v — > 0. 

Lemma 3.11 Under assumptions iTO)J , <TO)) . < T07]) . fiTWj) . fOg ) a™d P^?U - I3J3\) the solution 
S v to \3.27\) given by theorem \3.1U\ satisfies 



\VS u \' 2 dxd9) 1/2 < 



G 



(3.82) 



Ihr 



Proof of lemma [3.111 We proceed in a way similar to the proof of lemma \TM Multiplying (|3.27[) 
by S v and integrating in x over T 2 and in 6 over [0, 1] gives 



(A e + v)\\JS v \ 2 dxd6 < I I \C e \ \VS v \dxd0 < 7 



JT 2 



JT 2 



JT : 



the last inequality being obtained from (|3.13p . From (|3.83p we deduce 



'A. \VS" 



< 



7- 



L 2 (R,L 2 (T 2 )) 



'A IV^Idxdfl, (3.83) 



(3.84) 



On another hand since p.!2[) is assumed, we have 

7 IV7 C |2, 



G 



i/ir 



\ 1/2 

|V5 v | 2 da;d0 I 



T 2 



Sc. JT 2 



\ 1/2 






^/a|V5"| 





Formula (j3"^4l) and (j3~^5)) give (|3~g2"|) . 

As a direct consequence of lemma 13.111 we have the following corollary 



L|(R,L 2 (T 2 )) 

(3.85) 



COROLLARY 3.12 There exists 8 Q e [9 a ,6 u ] such that, under assumptions \3.8\) . \3.9\) . h3.ll}) , 
REUS), and [3~W - fOSp . S* satisfies 



|VS"(0o,OI|2 < 



(3.86) 



Gthr 



LEMMA 3.13 Under assumptions [3~8\) , l3~M . l3Al\) . [3~Wjl , l3~W\) and pJ^p - fUg) tfie solution 
S v of fr3.27\ l given by theorem \3.10\ satisfies 



IS* 



Ils : (r,l 2 (t 2 )) 



< 



•2 7 . 



(3.87) 
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Proof of lemma 13.131 First, because of p.75[) . doing the same as in the proof of lemma I3~5I we 
get from (|3.86l) that their exists a 9q such that: 

||$"(0o,-)ll2< -7==- (3-88) 

Gthr 

Secondly, in any 9 and any x where C e (9,x) ^ and A e (9,x) ^ 0, applying formula ()3.53j) with 
U = |V«S V | and V = \C e \ we obtain 

J ~ 2 



4 A e + v 

<^^|VSH 2 +7, (3.89) 

thanks to assumption p. 131) . Hence multiplying (|3.27p by S v and integrating over T 2 yields 

ld||S"0V)ll2 



dO 



(A + v)\VS"{9,-)\ 2 dx 

igT 2 , C t (6,x)=0 



+ I {Ae + v)\VS u (9,-)\ 2 dx 

JxeT 2 , C e (6,x)^0 

< [ \C e (0,x) -S7S u (9,-)\dx 

JxeT 2 , c c (0,x)^o 



/xeT 2 , C e (8,x)=£0 

< f ^^\VS»{9 1 -)\ 2 dx + / 1 dx. (3.90) 



Passing the first term of the right hand side in the left hand side yields 



dff 

Coupling (|3.88p and (|3.91[) allows us to deduce 



<2 7 . (3.91) 



S v (6,-)\\l < ~^= + ^l, (3.92) 



\J Gthr 

for any 9 € [9q,9q + 1] and because of the periodicity of S v , inequality (I3.87[) . 



LEMMA 3.14 Under assumptions [3~8\) . fOj) . ( TOT]) . i TOg)) . i TOffj) and pJ^| ) - [3~23\) the solution 
S v to l{3.27\ ) given by theorem l 3.10\ satisfies 



dt 



< 



L2f(R,L 2 (T 2 )) 



•27, 



G 



Ihr 



dS v 



< e 



,L 2 (T 2 )) 




(3.93) 



Proof of lemma 13.141 In a first step, remembering inequality (|3.84[) proved in the beginning of 
the proof of lemma 13.111 and using (|3.18|) , we deduce 



dt 



vs 1 - 



< 7 



,L 2 (T 2 )) 



i AXIS" 



<7 2 - 



(3.94) 



L 2 (R,L 2 (T 2 )) 
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In a second step, following the way to prove lemma |3~TI we multiply (|3.27p by — AS U and we integrate 
the resulting equality in x <E T 2 to get 



2 d6 



T 2 



{A e + v)\AS v \ 2 dx 



VA t ■ VS" AS u dx - V • C e AS"dx, (3.95) 
Jt 2 



which using (|3. 541) and (13.55)) and integrating in 6* over [0, 1] gives 



A e + V, 



lo Jt 2 2 
Now using (|3.16|) . we deduce 



1 f |VA|2 'VSTefcdfl 



o Jt 2 A e + v 



|VA| 2 
o Jt 2 A t + v 



o Jt 2 



Az\AS v \ z dxdO < / / A|VSTciEd6i + 7<7 2 +7 



o Jt 2 



cfcdtf, (3.96) 



(3.97) 



thanks to inequality (|3.84l) . Finally, because of (|3.17l) . we get from (|3.97[) 



at 



AS 1 



<7 



L|(RxT 2 ) 



'A, AS" 



L 2 (RxT 2 ) 



< 7 + Ta/7- 



(3.98) 



In the third step, we set out the equation to which \ s a solution. In a way similar to the one 



followed in the proof of lemma 13.91 we obtain 

/as* 

o 



where 



dt I f ~ dS u 

be -v-f(A(^r,,.) + ,)v— ) =v c 



dc € d_A 



(3.99) 



(3.100) 



Multiplying equation (13.99)) by and integrating in x G T 2 , in the same spirit as in the proofs of 



lemma I3T61 and \3M we deduce 
dS u 2 



8 



1 



2 90 



0f 



+ / (A + v) 

Jt 2 



,dS v 



at 



dx < 



T 2 



dC e 



at 



as v 



J2 



OA, 



dt 



at 
WS" 



dx 



7 dS u 
~~dt 



dx. 



To estimate the right hand side of (|3.101|) . we first notice that, applying ()3.23j) . we have 



/ 

Jt 2 



dC t 



dt 



,dS v 



at 



dx < 7 / \JA t 

T 2 



dS v 



at 



dx < 7 



dS v 



dt 



(3.101) 



(3.102) 



Then using p,17[) we deduce 





dA e 


\VS V \ 


as v 


dx < 




dAe 


/ 






It 2 


dt 


v at 






dt 



\vs v \ 



dAe 


as v 


at 


v at 



< 



A e NS v \ 



2 

dS v 
~dt 



(3.103) 
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As a consequence of (|3.102j) , (|3.103p and (|3~54"j) . integrating (|3.10ip in 9 over [0, 1] yields 

2 

< 7 VA 

L^^.L^T 2 )) L%(BL,L 2 (T 2 )) 



(A e + V) 



,dS u 
~dT 



7 dS u 



,dS v 
dt 



(3.104) 



L^(R,L 2 (T 2 )) 



where once again (|3.84p is also used. 
From this last inequality, we deduce 



,dS v 
dt 



L 2 (R,L 2 (T 2 )) 



< 7 + 7 



(3.105) 



and then 



,dS v 
dt 



7 + 7 



d9 < 
2 <'G 



(3.106) 



Ihr 



The fourth step consists in deducing from (|3.106p that there exists a 9q & [9 a ,9 u ] such that 



7 + 7" 



< 

2 'IG thT 



(3.107) 



and, since the mean value of ^§r{@0i ') is zero 



dS u 



dt 



'or, 



7 + 7 



< 

2 >'G 



(3.108) 



i/lr 



The fifth and last step, consists in going back to (|3.10ip . Applying formula (|3.53p with U = 
and V — |V^|-| to treat the first term of the right hand side of p.lOip and with U = \^\\VS U 



and V = |V^-| to treat the second, we get 



dS 1 * 



dt 



09 



T 2 2 



< 



,dS v 



dt 



dx 



dC e 


2 




dt 


—dx + / - 


dt 



IV5 



v I 2 



T 2 A t + V 



A e + V 



dx<j + A e \VS v \ z dx, (3.109) 



where we used (|3.23p and (|3.17p to find the last inequality. Integrating this last formula in 9 over 
[#0j cr ] for any a > 9q, we obtain, always remembering (|3.84p . 



dS u 



dt 



(<v) 



< 



dS "(9 ) 



o / 7 + 7 J 
2 */G 



+ 2 7 . 



(3.110) 



Inequality (|3. 1 10|) yields directly the first inequality of (|3.93p . using the periodicity of S". The proof 
of the second inequality of (|3.93p is done in a similar way. This ends the proof of lemma 13.141 ■ 
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As neither estimate (|3.87[) nor estimate (13.93[) depend on v, we can deduce that, extracting a subse- 
quence, as v — > 0, S v — > S in Lg(M., L 2 (T 2 )) weak-*, and that the limit satisfies estimates looking 
like (|3.87|) and ()3.93j) . but also a property of the type (|3.75|1 . and that it is solution to equation (|3.27|> 
with v — 0. In other words, we deduce the following theorem. 

Theorem 3.15 Under assumptions (EM > (SUM ) ULJM and WW - WIM) there exists a 

unique function S — S(t,T,0,x) G L°£(R, L 2 (T 2 )), periodic of period 1 with respect to 9, solution to 

dS ~ ~ 

— - V ■ (A e (t, T, ; -)V5) - V • C e (t, r, •, •), (3.111) 
and satisfying, for any t,T,9 € R + xlxl 

S(f, r, 0, i)cfe = 0. (3.112) 



Moreover it satsfies: 



T 2 



||^IU-( R> ^( T2 )) < -^== + 2 7 , (3.113) 



ds ,i2 ^ 7 + 7 



3 



II ^lli-(R,i»CP)) < ^7= + 2 7, (3.H4) 

V Gthr 

dS 112 / ^7 + 7 



ll^rlli-(K,^(T2)) < e rz + 2 7 y (3.115) 

Remark 3.1 Uniqueness of S is not gotten via the above evoked process v — > 0, but directly 
comes from (|3.111j) . Assuming that there are two solutions Si and 52 to (|3.111|) . we easily deduce 
that 

d(\\S 1 -S 2 \\l) , _ 

'-- \ A e \V(S 1 -S 2 )\ 2 dx = 0, (3.116) 

which gives, because of the non- negativity of A e , 



d (\\Si - S 



2 II 2 



d9 

From p,116p we deduce that either 



< 0. (3.117) 



or, for any 9 G 



A|V(5 1 -5 2 )| 2 = 0, (3.118) 

\\S X {9 + 1,.) - S 2 {9 + 1, 0||5 < \\Sx{d, ■) - S 2 {9, . (3.119) 

As p,119p is not possible because of the periodicity of Si and S 2 , we deduce that (|3.118[) is true. 
Using this last information, we deduce, for instance 

V(S 1 -S 2 )(9 u ,-) = 0, (3.120) 

yielding, because of property (|3.112p 

||(5! - S 2 )(0 U , OH" < l|V(Si - S 2 )(9 U , . (3.121) 
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Injecting (|3.118l) in (13.116)) yields 

d(\\Si-S»\\l) 

— =0, (3.122) 



de 

and then 

-5a)(«, Oils = 0, (3.123) 
for any 9 > 6 U and consequently or any 9e 1. ■ 
Theorem 1 1.1 1 is the consequence of theorem 13. 151 We will now end its proof. 

End of the proof of theorem [TTT1 Since, with definitions ([3TT ]) -(|3^g |l or (|3TT]) , ([3~3")) , (|3Td| . (f3~7|) . 
properties (JUIJ), $£9§, (|3~T2l) . (|3~T3)) and (j3~T4| - (j3~23f are consequences of assumptions ([L2]), (fl~3)) 
and (|1.4I) or p. 21) . (11.61) . (|1.7I) and (jl.8|) . theorem H~T1 is a direct consequence of the following theorem. 



Theorem 3.16 Under properties £0 )) . iT07)) . fOg) , ifQ3)) and pj^] ) - fgj^ ), for any T, 

not depending on e, equation i3.5]) . with coefficients given by \3.1\) (coupled with A3.2\) or \3.b}) ) 
and \3. 3\) (coupled with Jff.^| ) or \3. 7\j ), has a unique solution z e € L°°([0, T}; i 2 (T 2 )). This solution 
satisfies: 

lklz°°([ ,r],£ 2 (T 2 )) <7 (3.124) 
where 7 is a constant which does not depend on e. 

Proof of theorem 13.161 To prove uniqueness, we consider z\ and z\ two solutions of (|3.5|) . Their 
difference is then solution to 

9* e V^ V(Zl Z2 V~ U ' (3.125) 

( z l ~ z 2)|t=0 = ^' 

and multiplying the first equation of (I3.125[) by [z\ — zf) and integrating with respect to a: gives 

d(\\zl-zj\\l) 



dt 

yielding 



< 0, (3.126) 



||^-2!||=0, for any t, (3.127) 

and giving uniqueness. 

Concerning existence, corollary |3.2l gives existence of z e on a time interval of length of some e. 
Now we consider the function Z e = Z e (t, x) = S(t, -^=,^,x) where S is given by theorem 13. 151 Since 

3Z e dS . t t . 1 dS , t t . IdS. t t . 
■^7 = ^(t,-^,-,x) + -=-^-(t t -=,-,x) + --^(t,-=,-,x), 3.128 
at at ve e y/e or y/e e e 06 y/e e 

it is obvious to deduce from (|3.111|) that Z e is solution to 

dZ e 1 , Tr ,„„ 1„ , f 35, t i , IAS, it, ,„ „„, 

— V • AVZ E = -V • C £ + — (i, - 7=l + -p-^-Ct, -, x • 3.129 

at e e ot y/e e y/e or y/e e 

From (I3.129[) and (|3.5[) we deduce that z e — Z e is solution to 



d(z e -Z e ) 1 ,~ , f ^ \ \ 95, f i . 1 dS , t t . 
ot e at y/e e y/e or y/e e 

(z e -Z% =o = z o -S(0,0,-). 



(3.130) 



4 Homogenization, proof of theorem \1.2l 



25 



Multiplying (|3.130p by z e — Z e and integrating in x, using estimates (I3.114p and (|3.115[) gives: 



d(\\z<-Z<\\*) 
dt 



< 2 



~ + "' 3 + 2^\\z e -Z%, (3.131) 



\ \j Gthr 



which gives 



\\z e (t,-)-Z c (t,-)\\ 2 < 211^-5(0,0,0112 
which gives (|3.124[) with 



7 + 7' 



3 



•27*, (3.132) 



7 = 211*0-5(0,0,0112 ^^= + 2jT. (3.133) 

\ V Gthr 

Coupling local existence and estimate (|3. 124|) yield global existence and then theorem 13.161 is true. 

■ 

As a consequence theorem 1 1.1 1 is also true. ■ 

As a by-product of theorem l3.10l using a way similar of the one used to prove theorem 13. 161 we can 
obtain a theorem giving long-term existence of space-periodic solution to parabolic equation. 

Theorem 3.17 Under assumptions \3. 8]) . iS. 9\) and \3.11\) . for any v > 0, any e > 0, and any T not 
depending on e, equation \3.2J$ has a unique solution solution z e,l/ G L°°([0, T], L 2 (T 2 )). Moreover 

lk £ '1lL°=([o,T],L 2 (T 2 )) is bounded. 

4 Homogenization, proof of theorem 11.21 

We consider equation (|3.5p where A e and C e are defined by formulas (|3.1[) coupled with (|3.6p and 
(|3.3p coupled with (13. 7[) . Our aim consists in deducing the equations satisfied by the limit of z e 
solution to (13.51) as e — > 0. 



It is obvious that 

A e (t,x) two scale converges to A(t,0,x) € £°°([0, T], L# (R, L 2 (T 2 ))) 

and C c (i, x) two scale converges to C(t, 0, x), (4-1) 

with 

A(t, 9, x) = ag a (\U(t, 6,x)\) and C(t, d, x) = cg c (\U(t, 6, x)\) ffig *L (4.2) 
and we have the following theorem. 

Theorem 4.1 Under assumptions f£3p . nfflji . RTWj) . (MM) and W^W ~ WM, f or an V T , not 
depending on e, the sequence (z e ) of solutions to \3.5}) . with coefficients given by iS.l]) coupled with 
iTO|) and gnfy coupled with (3.7), two-scale converges to the profile U G L°°([0, T], L°£(R, L 2 (T 2 ))) 
solution to 

^L-V-(A\U)=V-C, (4.3) 



where A and C are given by J^.^p . 
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Proof of theorem 11.21 Theorem 11.21 is a direct consequence of theorem 14. II ■ 

Proof of theorem 14.11 Defining test function ip e (t,x) — tp(t,-,x) for any ip(t,6,x), regular 
with compact support in [0, T) x T 2 and periodic in 9 with period 1, multiplying (|3.5I) by -0 € and 
integrating in [0, T) x T 2 gives 

/ / ^r-^dtdx-- { { \7-(A e X7z e )iP e dtdx = -[ [ V-C e tp e dtdx. (4.4) 

Then integrating by parts in the first integral over [0,T) and using the Green formula in T 2 in the 
second integral we have 

zo(x)ip(0, 0, x)dx — [ [ —^—z e dtdx 
Jt 2 Jo 9t 

+ - I [ A e Vz e Vtp e dtdx = - / [ V ■ C c i) e dtdx. (4.5) 
6 Jt 2 Jo e Jt 2 Jo 

Again using the green formula in the third integral we obtain 



z o (x)ijj(0, 0, x) dx — I I —^j-z e dtdx 
Jt 2 Jo 

T 



at 



-II z e V ■ (A'W) dtdx = - I I V -C e ip e dtdx. 
e Jt 2 Jo e Jt 2 Jo 



But 



where 



dip 6 _ (dipY 1 (dip 
~dt ~ (~dt) + 1 ( ~m 



(4.6) 
(4.7) 



dipY , , dip, t . , (dipY , . dtp, t , 



\*((^)\ 1 -(^-Y+ 1 -v.(A*vr))dxdt 



t 2 Jo V \ dt J e \ 89 J e 

! , rT 



then we have 



V ■C t ip t dtdx = z (x)ip(0,Q,x)dx. (4.9) 
e Jt 2 Jo Jt 2 

Using the two-scale convergence due to Nguetseng [27J and Allaire [T] (see also Frenod Raviart and 
Sonnendriicker [13]), if a sequence / e is bounded in L°°(0, T, i 2 (T 2 )), then there exists a profile 
U(t,0,x), periodic of period 1 with respect to 9, such that for all ip(t,6,x), regular with compact 
support with respect to (t, x) and periodic of period 1 with respect to 9, we have 

PYdtdx — > III Utpd9dtdx, (4.10) 



t 2 Jo Jt 2 Jo Jo 



for a subsequence extracted from (f e ). 

Multiplying (|4.9[) by e and passing to the limit as e — > and using (|4.10p we have 



f I I U^d9dtdx + \im I / z e V ■ (A e V ip e ) dtdx 

T 2 J Jo dd c->°Jt 2 Jo 



lim / I C e ■ VYdtdx, (4.11) 
£^°Jt 2 Jo 
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for an extracted subsequence. As A e and C e are bounded (see (|3. 810 and ^ is a regular function, 
A e Vip e and \7ip e can be considered as test functions. Using (|4.ip we have 



/ / z £ V • (A e Vip e )dtdx — >[ [ [ UV ■ (AVi>)d9dtdx, (4.12) 
Jt 2 Jo Jt 2 Jo Jo 



and 

r T r rT r l 



[ [ C f - ■ Vip e dtdx two scale converges to f [ [ C ■ Vip dOdtdx. (4.13) 

Jt 2 Jo Jt 2 Jo Jo 

From this we obtain from (|4.11|) the equation satisfied by U : 

dU 



;m V-(AVU) =V-C, (4.14) 



which is equation (|4.3[) . 



Existence and uniqueness of equation (|4.3|) is given by theorem 13.151 (applied with e-independent 
coefficient and right hand side). 

From this uniqueness, we can deduce that the whole sequence (z e ) converges. 

Let us characterize the homogenized equation for A and C. We recall that A € and C e are given by 
formulas (|3.2[) and (|3.4[) . Multiplying this equation by a test function ip e using the two-scale limits 

A e ip e dtdx= [ [ a(l-6eA^(t,0,a;)) 3a (|W(t,6',a;)|)V e dida; 

T 2 Jo Jt 2 Jo 

ag a (\U(t,9,x)\)tlj e dtdx+ [ [ abeM{t,6,x) g a (\U(t,6,x)\)tp e dtdx. (4.15) 
It 2 Jo Jt 2 Jo 

Using the two-scale limits we have 

/ Aipd9dtdx= [ [ [ ag a (\U{t 7 9,x)\)i}dtdx. (4.16) 
It 2 Jo Jo Jt 2 Jo Jo 

Using the expression of C c we prove either that 



[ Ci/j dOdtdx = [ [ [ 
t 2 Jo Jo Jt 2 Jo Jo 



Crdfidhl.r = j j I rg c (\U(t,9,x)\) ^ ^ *j ^ dtdx, (4.17) 



then 

A = ag a (\U(t,9,x)\) and C = cg c (\U(t, 6, x)\) !f}!? i' % ■ (4.18) 

\u{t,e,x)\ 
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Considering (|3.5[) with coefficients given by (|3.ip coupled with (I3.6[) and p.3l) coupled with p.7p 
leads to writing 

A € (t,x) = A € (t,x) +eA{(t,x), (5.1) 

and 

C e (t,x) =C e {t,x) +eC{{t,x), (5.2) 
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where 

A e (t,x) = A(t, -,x), C e (t,x) =C(t, - x), (5.3) 
with A and C given by (|4.2I) and where 

Al(t,x)=Ai(t,-,x), C{(t,x)=C 1 (t,~,x), (5.4) 

with 

U(t,9,x) 
(5.5) 



Ai(t,9,x) = -abM(t,e,x) g a (\U(t,e,x)\) and d{t,e,x) = -cbM(t,9,x) g c (\U(t, 9, x)\) 

|W(t,0,a;)| 



Because of assumptions (|1.2p . (jl.3p and (jl.4|) . 

„4, C, Ai, Ci, A e , A\, C £ , and C\ are regular and bounded. (5.6) 
Moreover, the supplementary assumption that Uthr = yields 

A(t, 9, x) > Gthr for any t,0,x G [0, T]xRx T 2 . (5.7) 



Theorem 5.1 Under assumptions i3.8\) . i3.11\) . 115. 6\) and \5. 7| ), considering z e the solution to VS. 5\) 
with coefficients given by \3.1\) and i3. 3\) coupled with i5.1\) and \5. 2\) and U e = U e (t,x) = U(t,-,x) 
where U is the solution to \1.16\) . for any T not depending on e, the following estimate holds for 
z e - U e 

z e — U e 

— < a, (5.8) 

e L°°([0,T),L 2 (T 2 )) 

where a is a constant not depending on e. 

Furthermore, sequence ( z '~ u ' ) two-scale converges to a profile U 1 £ L°°([0, T], L°£(R, L 2 (T 2 ))) 
which is the unique solution to 



ATI 1 / ~ \ ~ FITJ ~ 

— -v\avu')=v-c 1 + — + v-(a 1 vu). 

Proof of theorem 11.31 Theorem ll.3l is a direct consequence of theorem 15. II 

Proof of theorem 15.11 Using (|5.1j) and (|5.2j) equations (|3.5j) becomes 

Bz e 1 ~ 1 ~ ~ ~ 

: V • {A c S7z e ) = -V • C e + V • (A\ Vz e ) + V • C\. 



dt 



For U e , we have 



where 



8U e 
~dt 



dU 
~dt 



1 fdU 



dU 
~dt 



(t,x) 



dU , t . 
— (t,-,x) and 
at e 



dU 
~d9 



, , dU , t . 



Using (|4.3p . U e is solution to 



(5.9) 



(5.10) 
(5.11) 

(5.12) 
(5.13) 
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Formulas (j5~TUl) and (|5TT3|) give 



dt e \ y 'J 1 \ at 

Multiplying equation (|5.14|l by - and using the fact that z c — z e — U e + U e in the right hand side 
of equation (15.141) . - ~ u is solution to: 



d 



-V .((A e + eAt)V(—--) = - V • C{ + ( — y + V • (A{VU<) . (5.15) 



dt e V v u v e V e 1 v dt 



Remark 5.1 Concerning notations, we have to pay attention to the fact that 

iVAandC e ^C £ . (5.16) 



Our aim here is to prove that z e u is bounded by a constant a not depending on e. For this let 
us use that _4 £ , A\ , C c and are regular and bounded coefficients (see (|5 .6|) ) and that .4 e > G t hr 
(see (J5J])). Thus, V • C{ is bounded, V • (A\VU' i ) is also bounded. Since U e is solution to (pTT3l) . 
^ satisfies the following equation 

Equation (|5.17l) is linear with regular and bounded coefficients. Then using a result of Ladyzenskaja, 
Solonnikov and Ural'Ceva [3T], ^ is regular and bounded . Then the coefficients of equations (I5.15[) 
are regular and bounded. Then, using the same arguments as in the proof of theorem 1 1.1 1 we obtain 
that (^f 1 ) is bounded, that it two-scale converges to a profile U 1 G L°°([0, T], L°g(R, i 2 (T 2 ))) 
and that this profile U 1 satisfies equation (|1.2ip . ■ 
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